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Abstract — We consider a dense, ad hoc wireless network, 
confined to a small region. The wireless network is operated 
as a single cell, i.e., only one successful transmission is 
supported at a time. Data packets are sent between source- 
destination pairs by multihop relaying. We assume that nodes 
self-organise into a multihop network such that all hops are 
of length d meters, where d is a design parameter. There is a 
contention based multiaccess scheme, and it is assumed that 
every node always has data to send, either originated from it 
or a transit packet (saturation assumption). In this scenario, 
we seek to maximize a measure of the transport capacity of 
the network (measured in bit-meters per second) over power 
controls (in a fading environment) and over the hop distance 
d, subject to an average power constraint. 

We first argue that for a dense collection of nodes confined 
to a small region, single cell operation is efficient for single 
user decoding transceivers. Then, operating the dense ad 
hoc wireless network (described above) as a single cell, we 
study the hop length and power control that maximizes the 
transport capacity for a given network power constraint. 
More specifically, for a fading channel and for a fixed 
transmission time strategy (akin to the IEEE 802.11 TXOP), 
we find that there exists an intrinsic aggregate bit rate (Q op t 
bits per second, depending on the contention mechanism and 
the channel fading characteristics) carried by the network, 
when operating at the optimal hop length and power control. 
The optimal transport capacity is of the form d op t(Pt) x & op t 

with d opt scaling as P t ~ , where Pt is the available time 
average transmit power and r\ is the path loss exponent. 
Under certain conditions on the fading distribution, we then 
provide a simple characterisation of the optimal operating 
point. 

Index Terms — Multihop Relaying, Self-Organisation 

I. Introduction 

We consider a scenario in which there is a large number 
of stationary nodes (e.g., hundreds of nodes) confined to 
a small area (e.g., spatial diameter 30m), and organised 
into a multihop ad hoc wireless network. We assume that, 
traffic in the network is homogeneous and data packets are 
sent between source-destination pairs by multihop relaying 
with single user decoding and forwarding of packets, i.e., 
signals received from nodes other than the intended trans- 
mitter are treated as interference. A distributed multiaccess 
contention scheme is used in order to schedule transmis- 
sions; for example, the CSMA/CA based distributed coor- 
dination function (DCF) of the IEEE 802.11 standard for 



wireless local area networks (WLANs). We assume that all 
nodes can decode all the contention control transmissions 
(i.e., there are no hidden nodes), and only one successful 
transmission takes place at any time in the network. In 
this sense we say that we are dealing with a single cell 
scenario. We further assume that, during the exchange of 
contention control packets, pairs of communicating nodes 
are able to estimate the channel fade between them and 
are thus able to perform power control per transmission. 

There is a natural tradeoff between using high power and 
long hop lengths (single hop direct transmission between 
the source-destination pair), versus using low power and 
shorter hop lengths (multihop communication using inter- 
mediate nodes), with the latter necessitating more packets 
to be transported in the network. The objective of the 
present paper is to study optimal routing, in terms of 
the hop length, and optimal power control for a fading 
channel, when the network (described above) is used in 
a multihop mode. Our objective is to maximise a certain 
measure of network transport capacity (measured in bit- 
meters per second; see Section [TVT >. subject to a network 
power constraint. A network power constraint determines, 
to a first order, the lifetime of the network. 

Situations and considerations such as those that we 
study could arise in a dense ad hoc wireless sensor 
network. Ad hoc wireless sensor networks are now being 
studied as possible replacements for wired measurement 
networks in large factories. For example, a distillation col- 
umn in a chemical plant could be equipped with pressure 
and temperature sensors and valve actuators. The sensors 
monitor the system and communicate the pressure and 
temperature values to a central controller which in turn 
actuates the valves to operate the column at the desired op- 
erating point. Direct communication between the sensors 
and actuators is also a possibility. Such installations could 
involve hundreds of devices, organised into a single cell ad 
hoc wireless network because of the physical proximity of 
the nodes. There would be many flows within the network 
and there would be multihopping. We wish to address the 
question of optimal organisation of such an ad hoc network 
so as to maximise its transport capacity subject to a power 
constraint. The power constraint relates to the network life- 



time and would depend on the application. In a factory 
situation, it is possible that power could be supplied to 
the devices, hence large power would be available. In 
certain emergencies, "transient" sensor networks could 
be deployed for situation management; we use the term 
"transient" as these networks are supposed to exist for 
only several minutes or hours, and the devices could be 
disposable. Such networks need to have large throughputs, 
but, being transient networks, the power constraint could 
again be loose. On the other hand sensor networks de- 
ployed for monitoring some phenomenon in a remote area 
would have to work with very small amounts of power, 
while sacrificing transport capacity. Our formulation aims 
at providing insights into optimal network operation in a 
range of such scenarios. 

A. Preview of Contributions 

We motivate the definition of the transport capacity of 
the network as the product of the aggregate throughput 
(in bits per second) and the hop distance (in meters). For 
random spatio-temporal fading, we seek the power control 
and the hop distance that jointly maximizes the transport 
capacity, subject to a network average power constraint. 
For a fixed data transmission time strategy (discussed in 
Section UlI-BI) . we show that the optimal power allocation 
function has a water pouring form (Section IV-AI ). At the 
optimal operating point (hop distance and power control) 
the network throughput (O op t, in bits per second) is 
shown to be a fixed quantity, depending only on the 
contention mechanism and fading model, but independent 
of the network power constraint (Section IV-Bb . Further, 
we show that the optimal transport capacity is of the form 

dopt(Pt) x ©opt, with d opt scaling as Pt n , where P t is the 
available time average transmit power, and r\ is the power 
law path loss exponent (Theorem l5.2b . Finally, we provide 
a condition on the fading density that leads to a simple 
characterisation of the optimal hop distance (Section lV-Ct . 

II. Motivation for Single Cell Operation 

In this context (a dense, ad hoc wireless network), the 
seminal paper by Gupta and Kumar |]3] would suggest that 
each node should communicate with neighbours as close 
as possible while maintaining network connectivity. This 
maximises network transport capacity (in bit-meters per 
second), while minimising network average power. It has 
been observed by Dousse and Thiran in [4|, that if, unlike 
0, a practical model of bounded received power for finite 
transmitter power is used, then the increasing interference 
with an increasing density of simultaneous transmitters is 
not consistent with a minimum SINR requirement at each 



receiver. The following discussion motivates single cell 
operation for our framework. 

Consider a planar wireless network with n nodes in a 
square of fixed area A. Let K(n) be the spatial reuse in the 
network (the number of simultaneous transmissions) and 
r(K(n)), the maximum transmitter-receiver separation. 
Denote by P(K(n)) the common power per transmitter 
(assumed to be the same for all nodes) satisfying a 
network average power constraint P (as in Section fVt . 
The maximum SINR achievable per link in such a network 
(with single user decoding receivers) is bounded above 
by (jv+jj^i ' wnere N is the receiver noise power and 



l K(n): 



the interference at a node due to spatial reuse. 



Using the finite (and fixed) area assumption, the mini- 
mum interference power from any simultaneous transmis- 
sion is given by P W"^ Hence, the SINR achievable 

P(K(n)) 



P(K(n)) 

over any link is bounded above by 



(2A) 2 

and the maximum bit rate achievable over a link is, 



log 1 



P(K(n)) 
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| . The aggregate throughput 



N+(K(n)- 

in the network is now bounded above by C(K(n)), 



C(K(n)) :=K(n)log 1 + 



P(K(n)) 



N+(K(n) - 1) 



P(K(n)) 



Clearly, C(K(n)) is uniformly bounded aboveQ. Also, 
r(K(n)) < \/2A. Hence, we see that the transport 
capacity achievable in the network (bounded above by 
^V{K{n):K(n)>2}C{K{n))r{K{n))) is finite, indepen- 
dent of the number of nodes and the network power P. 
Further, we would expect the transmitter-receiver sepa- 
ration (bounded above by r(K(n))) to decrease to as 
K(n) increases to oo (finite area assumption). Hence, 
limj ! -(„)^ 00 C(K(n))r(K(n)) = 0. This implies that 
there exists an optimal K(n),l < K(n) < oo, which 
maximises the transport capacity in the network, i.e., 
the optimum spatial reuse is finite. Now, consider a 
simple TDMA scheme, without spatial reuse and with 
direct transmission between the source and the destination. 
For a transmit power P, the TDMA schedule achieves 
0(log(P)) transport capacity, i.e., an unbounded transport 
capacity as a function of the network power P. As 
discussed above, with spatial reuse, the system becomes 
interference limited, and hence, becomes inefficient for 
large P. More recently, El Gamal and Mammen |5| have 
shown that, if the transceiver energy and communication 

'Note that C(K(n)) < K(n)\og ^1 + p^jxj Y independent of 
the transmit power P(K(n)). 



overheads at each hop is factored in, then the operating 
regime studied in Q is neither energy efficient nor delay 
optimal. Fewer hops between the transmitter and receiver 
(and hence, less spatial reuse) reduce the overhead energy 
consumption and lead to a better throughput-delay tradeoff. 

While optimal operation of the network might suggest 
using some spatial reuse (finite, as discussed above), 
coordinating simultaneous transmissions (in a distributed 
fashion), in a constrained area, is extremely difficult and 
the associated time, energy and synchronisation overheads 
have to be accounted for. In view of the above discussion, 
in this paper, we assume that the medium access control 
(MAC) is such that only one transmitter-receiver pair 
communicate at any time in the network. 

A. Outline of the Paper 

In Section Hn] we describe the system model and in Sec- 
tion |IV] we motivate the objective. We study the transport 
capacity of a single cell multihop wireless network, oper- 
ating in the fixed transmission time mode, in Section IV1 
Section PvTI concludes the paper and discusses future work. 

III. The Network Model 

There is a dense collection of immobile nodes that 
use multiaccess multihop radio communication with single 
user decoding and packet forwarding to transport packets 
between various source-destination pairs. 

• All nodes use the same contention mechanism with 
the same parameters (e.g., all nodes use IEEE 802.11 
DCF with the same back-off parameters). 

• We assume that nodes send control packets (such as 
RTS/CTS in IEEE 802.11) with a constant power 
(i.e., power control is not used for the control pack- 
ets) during contention, and these control packets 
are decodable by every node in the network. As 
in IEEE 802.11, this can be done by using a low 
rate, robust modulation scheme and by restricting the 
diameter of the network. This is the "single cell" 
assumption, also used in 0], and implies that there 
can be only one successful ongoing transmission at 
any time. 

• During the control packet exchange, each transmit- 
ter learns about the channel "gain" to its intended 
receiver, and decides upon the power level that is 
used to transmit its data packet. For example, in 
IEEE 802.11, the channel gain to the intended re- 
ceiver could be estimated during the RTS/CTS control 
packet exchange. Such channel information can then 
be used by the transmitter to do power control. In 
our paper, we assume that such channel estimation 



and power control is possible on a transmission-by- 
transmission basis. 

• In this work, we model only an average power con- 
straint and not a peak power constraint. 

• Saturation assumption : We assume that the traffic 
is homogeneous in the network and all the nodes 
have data to send at all times; these could be lo- 
cally generated packets or transit packets. In Q, 
the authors study the problem of load balancing 
in dense multihop wireless networks with arbitrary 
traffic requirements. In our work, we do not restrict 
to straight line paths, and permit such a load balancing 
routing strategy as in Q, which ensures that the load 
and the channel access pattern are identical for all the 
nodes. 

Data packets are sent between source-destination pairs by 
multihop relaying. Based on the dense network and traffic 
homogeneity assumption, we further make the following 
assumption. 

• The nodes self-organise so that all hops are of length 
d, i.e., a one hop transmission always traverses a 
distance of d meters. This hop distance, d, will be 
one of our optimisation variables. 

For a random node deployment, the hop distance that 
maximizes the system throughput need not be the same 
for every node and every flow. However, the approxima- 
tion holds good for a homogeneous network with large 
number of nodes. Further, it will be practically infeasible 
to optimize every hop in a dense setup with hundreds of 
nodes. 

A. Channel Model: Path Loss, Fading and Transmission 
Rate 

The channel gain between a transmitter-receiver pair 
for a hop is assumed to be a function of the hop length 
(d) and the multipath fading "gain" (h). The path loss 
for a hop distance d is given by -4-, where 77 is the 
path loss exponent, chosen depending on the propagation 
characteristics of the environment (see, for e.g., [8|). This 
variation of path loss with d holds for d > do, the far 
field reference distance; we will assume that this inequality 
holds (i.e., d > do), and will justify this assumption in the 
course of the analysis (see Theorem 15. 2\ . 

We assume a flat and slow fading channel with additive 
white Gaussian noise of power a 2 . We assume that for each 
transmitter-receiver pair, the channel gain due to multipath 
fading may change from transmission to transmission, but 
remains constant over any packet transmission duration. 
Since successive transmissions can take place between 
randomly selected pairs of nodes (as per the outcome 
of the distributed contention mechanism) we are actually 



modeling a spatio-temporal fading process. We assume that 
this fading process is stationary in space and time with 
some given marginal distribution H. Let the cumulative 
distribution of H be A(h) (with a p.d.f. a(h)), which by 
our assumption of spatio-temporal stationarity of fading 
is the same for all transmitter-receiver pairs and for all 
transmissions. We assume that the channel coherence time, 
t c , applicable to all the links in the network, upper bounds 
every data transmission duration in the network. Further, 
we assume that H and r c are independent of the hop 
distance d. 

When a node transmits to another node at a distance d 
(in the transmitting antenna's far field), using transmitter 
power P, with channel power gain due to fading, h, then 
we assume that the transmission rate given by Shannon's 
formula is achieved over the transmission burst; i.e., the 
transmission rate is given by 

where W is the signal bandwidth and a is a constant ac- 
counting for any fixed power gains between the transmitter 
and the receiver. Note that this requires that the transmitter 
has available several coding schemes of different rates, one 
of which is chosen for each channel state and power level. 

B. Fixed Transmission Time Strategy 

We consider a fixed transmission time scheme, where 
all data transmissions are of equal duration, T (< r c ) sees, 
independent of the bit rate achieved over the wireless link. 
This implies that the amount of data that a transmitter 
sends during a transmission opportunity is proportional to 
the achieved physical link rate. Upon a successful control 
packet exchange, the channel (between the transmitter, 
that "won" the contention, and its intended receiver) is 
reserved for a duration of T seconds independent of the 
channel state h. This is akin to the "TxOP" (transmission 
opportunity) mechanism in the IEEE 802.11 standard. 
Thus, when the power allocated during the channel state 
h is P(h), C(h)T bits are sent across the channel, 
where C{h) = Wlog (l + ^^). When P(h) = 0, 
we assume that the channel is left idle for the next T 
seconds. The transmitter does not relinquish the channel 
immediately, and the channel reserved for the transmitter- 
receiver pair (for example, by the RTS/CTS signalling) is 
left empty for the duration of T seconds. 

The optimality of a fixed transmission time scheme, for 
throughput, as compared to a fixed packet length scheme, 
can be formally established (see Appendix IST l. we only 
provide an intuition here. When using fixed packet lengths, 
a transmitter may be forced to send the entire packet even 



if the channel is poor, thus taking longer time and more 
power. On the other hand, in a fixed transmission time 
scheme, we send more data when the channel is good and 
limit our inefficiency when the channel is poor. 

IV. Multihop Transport Capacity 

Let d denote the common hop length and {P(h)} a 
power allocation policy, with P(h) denoting the transmit 
power used when the channel state is h. We take a simple 
model for the random access channel contention process. 
The channel goes through successive contention periods. 
Each period can be either an idle slot, or a collision period, 
or a successful transmission with probabilities Pi,p c and 
p s respectively. Under the node saturation assumption, the 
aggregate bit rate carried by the system, 0t({P(/i)}, d), 
for the hop distance d and power allocation {P(h)}, is 
given by (see U) 

p.(jz°m dA(h)) 



e T ({P(h)},d) := 



(1) 



PiTi + p c T c + p s (T + T) 

where L(h) := C(h)T (and C{h) is a function of {P(h)} 
and d), Ti,T c and T a are the average time overheads 
associated with an idle slot, collision and data transmis- 
sion. For example, in IEEE 802.11 with the RTS/CTS 
mechanism being used, a collision takes a fixed time 
independent of the data transmission rate. We note that 
Pi,PstPc,Ti,T , and T c depend only on the parameters of 
the distributed contention mechanism (MAC protocol) and 
the channel, and not on any of the decision variables that 
we consider. 

With 0-r({P(ft.)}, d) defined as in ((T), we consider 
0y({P(/i,)}, d) x d as our measure of transport capacity 
of the network. This measure can be motivated in several 
ways. 0y({P(/i)}, d) is the rate at which bits are transmit- 
ted by the network nodes. When transmitted successfully, 
each bit traverses a distance d. Hence, 0r({P(/i)}, d) x d 
is the rate of spatial progress of the flow of bits in the 
network (in bit-meters per second). Viewed alternatively, it 
is the weighted average of the end-to-end flow throughput 
with respect to the distance traversed. Suppose that a flow 
i covers a distance Di with ^ hops (assumed to be an in- 
teger for this argument). Let /3j©x({P(/i)}, d) be the frac- 
tion of throughput of the network that belongs to flow i. 
Then, Si^zi^lMhBl j s me end-to-end throughput for flow 

i and f^T({P(h)}4) xD . = / g i e T ({p(/ l )}, d) x d is the 

d 

end-to-end throughput for flow i in bit-meters per second. 
Summing over all the flows, we have 0x({P(/i)}, d) x d, 
the aggregate end-to-end flow throughput in bit-meters per 
second. 

With the above motivation, our aim in this paper is to 
maximise the quantity 9y({P(/i)}, d) x d over the hop 



distance d and over the power control {P(h)}, subject to 
a network average power constraint, P. We use a network 
power constraint that accounts for the energy used in 
data transmission as well as the energy overheads asso- 
ciated with communication. The network average power, 
V{{P{h)}), is given by, 

V({P(h)}) := 

Pi Ei + Pc E c + Ps (E + T J °° P(h) dA(h) ) 
PiTi +p c T c + p s (T + T) 
Ei , E c and E a correspond to the energy overheads as- 
sociated with an idle period, collision and successful 
transmission. Thus, Ei denotes the total energy expended 
in the network over an idle slot, E c denotes the total 
average energy expended by the colliding nodes, as well 
as the idle energy of the idle nodes, and E denotes the 
average energy expended in the successful contention ne- 
gotiation between the successful transmitter-receiver pair, 
the receive energy at the receiver (in the radio and in the 
packet processor), and the idle energy expended by all the 
other nodes over the time T a + T. We assume that Ei , E c 
and E D depend only on the contention mechanism and not 
on the decision variables d and {P(h)}. 

V. Optimising the Transport Capacity 

For a given {P(h)} and d, and the corresponding 
throughput <dT({P(h)},d), the transport capacity in bit- 
meters per second, which we will denote by -0({P(/i)}, d), 
is given by 

tP({P(h)},d) :=e T ({P(h)},d)xd 

Maximizing ip(-,-) involves optimizing over d, as well 
as {P(h)}. However, we observe that, it would not be 
possible to vary d with fading, as routes cannot vary at 
the fading time scale. Hence, we propose to optimize first 
over {P(h)} for a given d, and then optimize over d, i.e., 
we seek to solve the following problem, 

max max ifj({P(h)}, d) (3) 

d {{P(h)}:V({P(h)})<P} 

For a given d and power allocation {P(h)}, define the 
time average transmission power, P t ({P(h)},d), and the 
time average overhead power, P Q , as 

p s (J™P(h) dA(h) )T 



P t ({P(h)},d) := 



PiTi + p c T c + p s (T + T) 

p i E l +p c E c +p s E 
PiT + p c T c + p s (T + T) 



Observe that P a does not depend on {P(h)} and d. Now, 
the network power constraint can be viewed as 

Pt({P(h)},d)<P-P 



where the right hand side is independent of {P(h)} or 
d. Pt := P — P , is the time average transmitter power 
constraint for the network. 

A. Optimization over {P(h)} for a fixed d 
Consider the optimization problem (from (01) 

max ip({P(h)},d) (4) 

{{P(h)}:V({P(h)})<P} 

The denominators of 0t(-,-) in (fl]l and of V in (f2]i 
are independent of d and the power control {P(h)}. 
Thus, with d fixed, the optimization problem simplifies 
to maximizing J Q L(h) dA(h) or, 



log 



, P(h)ha\ , .... 

1 + ^r) dA{h) 



subject to the power constraint, 



P(h) dA(h) < Pt 



where P t is given by, 

pi , = (PtT% + PcT c + p s (T + T)) p 
p s T 

Pt is the average transmit power constraint averaged only 
over the transmission periods (successful contention slots). 

Without a peak power constraint, this is a well-known 
problem whose optimal solution has the water-pouring 
form (see O). The optimal power allocation function 
{P(h)} is given by 



P(h) 



1 dV 2 ^ + 



A ha 

where A is obtained from the power constraint equation 



a(h)P(h)dh = P t ' 



The optimal power allocation is a nonrandomized policy, 
where a node transmits with power P(h) every time the 
channel is in state h (whenever P(h) > 0), or leaves the 
channel idle for h such that P(h) = 0. 

B. Optimization over d 

By defining := the problem of maximising 

the throughput over power controls, for a fixed d, becomes 

max J log ( 1 H fCC 1 ) ) a(h)dh 

subject to 

/■OO p i 

t\(h)a(h)dh < JL 



Observe that P[ and d influence the optimization problem 
only as ||-. Denoting by T LprJ the optimal value of this 
problem, the problem of optimisation over the hop-length, 
d, now becomes 



max a x 1 — — 
d \ do 



(5) 



Theorem 5.1: In the problem defined by (0, the ob- 
jective d X r ("pr)' when viewed as a function of d, is 
continuously differentiable. Further, when the channel fad- 
ing random variable, H, has a finite mean (E(H) < oo), 
then 

1) lim d ^o d X T (f^) =0and, 

2) if in addition, rj > 2, r?d (r) is continuously 
differentiable and P(iZ > h) = O (p-) for large 

/i, then, lim^oo dxT (Sfpj = 0, 
Proof: The proofs of continuous differentiability of 
d x r r§j")> 1) an d 2) are provided in Appendix |B1 ■ 
Remarks 5.1: 

1) Under the conditions proposed in Theorem 15.11 it 
follows that d x V is bounded over d G [0, oo) 
and achieves its maximum in d G (0, oo). 

2) When the objective function (0 is unbounded, the 
optimal solution occurs at d — oo (follows from the 
continuity results). 

3) We note that, in practice, r\ > 2. 

Let do be the far field reference distance (discussed in 
Section IHLA1 . 

Theorem 5.2: The following hold for the problem in 

©. 

1) Without the constraint d > do, the optimum hop 
distance d opt scales as (P t )v. 

2) There is a value Pt min sucn that, for P t ' > Pt' m i n , 
d op t > do, and the optimal solution obeys the scaling 
shown in 1). 

3) For Pt > Amin> me optimum power control 
{P(h)} is of the water pouring form and scales as 
Pt- _, 

4) For Pt > Ptpu n , the optimal transport capacity 
scales as (iV)''. 

Proof: 

1) Let d op t be optimal for P\ > 0. We claim that, for 
x > 0, x~d op t is optimal for the power constraint 
xPt ■ For suppose this was not so, it would mean 
that there exists d > such that 



:i d, 



xR 



opt 



(X^dopt) 11 / 



<dT 



xP t 
dP 



or, equivalently, 

/ Pt' 



d pt r 



a opt 



< x 'd r 



Pt 



\(x~~d) T i J 

which contradicts the hypothesis that d op t is optimal 
for P t '. 

2) With the path loss model we see that for d < do, 
the received power is scaled more than the transmit- 
ted power P, due to the factor -4-, and an factor in 
a, i.e., the model over-estimates the received power 
and the transport capacity. Hence, the achievable 
transport capacity for d < do is definitely less than 
dx r {%rj ■ The result now follows from the scaling 
result in 1). 

3) It follows from 1) that, if Pt scales by a factor 
x, then the optimum d scales by x*i, so that, at 
the optimum, is unchanged. Hence the optimal 

is unchanged, which means that {P(h)} 
must scale by x. The water pouring form is evident. 

4) Again, by 1) and 2), if P t scales by a factor x, then 
the optimum d scales by x~ , so that, at the optimum, 
^j- is unchanged. Thus T (^^rj is unchanged, and 
the optimal transport capacity scales as the optimum 
d, i.e., by the factor xn. 

■ 

Remarks 5.2: 

The above theorem yields the following observations for 
the fixed transmission time model. 

1) As an illustration, with rj = 3, in order to double 
the transport capacity, we need to use 2 3 times the 
Pt . This would result in a considerable reduction in 
network lifetime, assuming the same battery energy. 

2) We observe that as the power constraint Pt scales, 
the optimal bit rate carried in the network, r ( Sfc- J , 
stays constant, but the optimal transport capacity 
increases since the optimal hop length increases. 
Further, because of the way the optimal power 
control and the optimal hop length scale together, 
the nodes transmit at the same physical bit rate in 
each fading state; see the proof of Theorem |5.2| part 
3). ■ 

C. Characterisation of the Optimal d 

By the results in Theorem 15.11 we can conclude that 
the optimal solution of the maximisation in (0 lies in the 
set of points for which the derivative of x T ^ j is 

zero. For a fixed P t ', define ir(d) := -p r . Differentiating 
d x r(7r(d)), we obtain, (see Appendix [Ab 

d 



dd 



(d r(7r(d)) 



r(7r(<Q) - r]n{d)X(ir(d)) 




Fig. 1, Plot of d X r y-jzj (lineal' scale) vs. d (log scale) for a channel 
with two fading states h\,h%. The fading gains are h\ = 100 and 
h>2 = 0.5, with probabilities a%. = 0.01 = 1 — a,h 2 - The function has 
3 non-trivial stationary points. 



where X(tt) is the Lagrange multiplier for the optimisation 
problem that yields r(7r(<i)). Since d appears only via 
7r(d), we can view the right hand side as a function 
only of 7r. We are interested in the zeros of the above 
expression. Clearly, it = is a solution. The solution 
7r = corresponds to the case d = oo; However, we are 
interested only in solutions of d in (0,oo), and hence, we 
seek positive solutions of it of 



r(-7r) — ry7rA(7r) = 



(6) 



Remarks 5.3: In Appendix [A] we consider a continu- 
ously distributed fading random variable H with p.d.f. 
a(h). The analysis can be done for a discrete valued 
fading distribution as well, and we provide this analysis in 
Appendix Icl The following example then illustrates that, 
in general, the function r(7r) — r)ir\(Tr) — can have 
multiple solutions. Consider a fading distribution that takes 
two values: hi = 100 and /12 = 0.5, with probabilities 
a/ ll = 0.01 = 1 — a/i 2 . Figure Q] plots d x T (^3-) for 
the system with r\ = 3. Notice that there are 3 stationary 
points other than the trivial solution d = 00 (which is not 
shown in the figure). Also, the maximising solution is not 
the first stationary point (the stationary point close to 0). 
If, on the other hand, = 0.001 = 1 — a,h 2 , we again 
have 3 stationary points, but the optimal solution now is 
the first stationary point. ■ 

More generally, and still pursuing the discrete case, let 
Ti denote the set of fading states when the fading random 
variable is discrete with a finite number of values; \H\ 
denotes the cardinality of 7i. 

Theorem 5.3: There are at most 2\H\ — 1 stationary 
points of d T(it(d)) in < d < 00. 

Proof: See Appendix [C] for the related analysis and 



the proof of this theorem. ■ 
We conclude from the above discussion that it is dif- 
ficult to characterise the optimal solution when there are 
multiple stationary points. Hence we seek conditions for 
a unique positive stationary point, which must then be the 
maximising solution. In Appendix [A] we have shown that 
the equation characterising the stationary points, r(-7r) — 
?77rA(7r) = 0, can be rewritten as 

[ (l og (y)-r,(y-l))^f(-)dy = (7) 

Jo y w j 



for f(x) :— a {^-^j 77 > tne density of the random 
variable X :— ^pf. Notice that n does not appear in 
this expression. The solution directly yields the Lagarange 
multiplier of the throughput maximisation problem for 
the optimal value of hop length. The following theorem 
guarantees the existence of at most one solution of (0. 

Theorem 5.4: If for any Ai > A 2 > ^ 



is 



a strictly monotonic decreasing function of y, then the 
objective function dxT (^^rj nas at most one stationary 
point d opU < d opt < 00. 

Proof: The proof follows from Lemmas IA.1I and IA.2I 
in Appendix lAl ■ 
Corollary 5.1: If H has an exponential distribution and 
77 > 2, then the objective in the optimisation problem of 
(O has a unique stationary point d op t E (0, 00), which 
achieves the maximum. 

Proof: a(h) is of the form /ie^' 1 . From Theorem l5.ll 
we see that lim^o d x T (^^r^j = an d lim^oo d x 
r f ^~J = 0. And, the monotonicity hypothesis in Theo- 
rem 15 .41 holds for a(h). ■ 
Remarks 5.4: 1) Hence, for 77 > 2, for the Rayleigh 
fading model there exists a unique stationary point 
which corresponds to the optimal operating point. 
2) For P t > -Ptmini an d for the conditions in Theo- 
rem [5J] and [54] let ir opt denote the unique stationary 
point of ©. Then define r(7T op t) = <d pt- It follows 
from Theorem [572] that the optimal transport capacity 



takes the form 



B op t, where O opt depends 



on a(h) and the MAC parameters but not on P (or 
Pt). 

3) Figure [2] numerically illustrates our results for the 
Rayleigh fading distribution and 77 = 2. Scaling P t 
by 4 scales the transport capacity from 2.3 to 4.6, 
i.e., by 4^ — V4 and similarly for scaling Pt by 9. 

The uniqueness result guarantees that a distributed im- 
plementation of the optimization problem, if it converges, 




10" 2 10" 1 10° 10' 10 2 



Fig. 2. Plot of d X r ( ^pfj (linear scale) vs. n (= -^-) (log scale) for 
a fading channel (with exponential distribution). We consider 3 power 
levels (Pt , 4Pt and 9Pt ) and rj = 2. The function has a unique 
optimum ir op t(iT op t ~ 0.2) for all the 3 cases. 

shall converge to the unique stationary point, which is the 
optimal solution. 

VI. Conclusion 

In this paper we have studied a problem of optimal 
power control and self-organisation in a single cell, dense, 
ad hoc multihop wireless network. The self-organisation is 
in terms of the hop distance used when relaying packets 
between source-destination pairs. 

We formulated the problem as one of maximising the 
transport capacity of the network subject to an average 
power constraint. We showed that, for a fixed transmission 
time scheme, there corresponds an intrinsic aggregate 
packet carrying capacity at which the network operates 
at the optimal operating point, independent of the aver- 
age power constraint. We also obtained the scaling law 
relating the optimal hop distance to the power constraint, 
and hence relating the optimal transport capacity to the 
power constraint (see Theorem 15.2b . Because of the way 
the power control and the optimal hop length scale, the 
optimal physical bit rate in each fading state is invariant 
with the power constraint. In Theorem 15.41 we provide a 
characterisation of the optimal hop distance for cases in 
which the fading density satisfies a certain monotonicity 
condition. 

One motivation for our work is the optimal operation 
of sensor networks. If a sensor network is supplied with 
external power, or if the network is not required to have 
a long life-time, then the value of the power constraint, 
P, can be large, and a long hop distance will be used, 
yielding a large transport capacity. On the other hand, if 
the sensor network runs on batteries and needs to have 
a long life-time then P would be small, yielding a small 



hop length. In either case, the optimal aggregate bit rate 
carried by the network would be the same. 

In (6), the author studies the problem of developing 
a distributed algorithm for nodes to adapt themselves 
towards the optimal operating point. They first propose a 
distance discretization technique in which the hop distance 
on the critical geometric graph is used as a distance 
measure. Using the distance approximation, the author 
then develops a distributed algorithm aimed to maximize 
the transport capacity of the network in the sense of the 
framework presented in this paper. 
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Appendix 

A. Stationary Points of d x T{i:{d)) 

Recall that we defined n(d) := Further, r(7r(d)) 
was defined as 

r(7r(d)) := max £ log + <h)dh (8) 

where the maximum is over all power controls {P(h)} 
satisfying the constraint 

J ^a(h)dh < ir(d) (9) 

For ease of notation, let us use the substitution x :— 
Write i{x) := = ^ and f(x) := a \. 

Note that /(•) is the probability density of the random 
variable X := Then, equations ([8]) and (O can be 
rewritten as 

/>oo 

r(7r) = max / log(l + x^{x))f{x)dx 
Jo 



and 



£(x)f(x)dx < it 



This optimisation problem is one of maximising a convex 
functional of {£(x)}, subject to a linear constraint. The 
optimal solution of the problem has water-pouring form, 
and the optimal solution is given by, 



i 



v A(7r) X 
where A(ir) is obtained from 



I f(x)dx = 7T 



a(tt) \A(tt) X 

Further, the derivative of the optimum value r(7r), w.r.t. 
7T, i.e., ^1 = A(tt) (see Aubin 0). 

Let us now reintroduce the dependence on d, and 
consider the problem of maximizing d x r(7r(d)) over d. 
Differentiating d x T(ir(d)) w.r.t. d, we get, 

^(dT(n(d))) = T(n(d))+d^T(n(d)) 

= r M + dfwd))xM 



= r(?r(d)) +rf r'(7r(d)) x 



-r?P f 



= r(7r(d)) - ij7r(d)r'(7r(d)) 

where r'(7r) := Substituting r(?r) = A(tt), we 

have, 

^(dT(7r(d))) = r(7r(d)) - »77r(d)A(7r(d)) (10) 

The stationary points of d x r(ir(d)) are now obtained 
by equating the right hand side of ( fTOb to zero. Note that 
since d appears in this equation only as n(d), we need 
only study the roots of the equation 



r(7r) — rjir\(ir) = 



(11) 



We now proceed to obtain a characterisation of the 
stationary points. Substituting the optimal solution in the 
expression of r(7r) and A(7r), and suppressing the argu- 
ment 7r in A(7r), we get, 

1» = J" log (I) f( X )dx (12) 

with A being given by 

* = [ (i - i) f{x)dx (13) 



Using the substitution z — ^, I — j, and defining g(^) 
?/(7),0and0 becomes, 



r(Tr) 



log 



g(z)dz 



(14) 



(15) 



/o 

with I (actually, being given by 
[* 

n = (Z — z) g(z)dz 
Jo 

We note that g(-) is the density of the random variable 

7 j_ _ _of_ 

^ -— X — off' 

We will use the following definitions for convenience. 
For a function t(-) of the random variable Z, define the 
operators E/(-) and G/(-) as 



So t{z)g{z)dz 
So g(z)dz 



Gi(t(Z)) := [ t(z)g(z)dz 
Jo 

Lemma A.l: The roots of ( fTTT i are equivalent to the 
roots of the equation 



I 



1 = G, log 



(16) 



with I then being given by ( TT~5T >. 

Proof: Using the definitions of E;(-) and G/(-), ( TBI 
and (TT~5T > simplify to 

r(7r) = log(l)P(Z < - G,(log(Z)) (17) 

7T = iP(Z < - Gi(Z) (18) 

( fT8l provides the I (actually 1(tt)) to be substituted in ( TT7| >. 
Substituting for r(7r) (from and for I (from (fT8l), 

into (fTTt . dividing across by P(Z < Z), and using the 
definition of E/(-), we have, 

'iT + Gtizy 



log 



P(Z < I) 



Ej(log(Z)) 



1]H 
7T + Gf(Z) 







log 



P(Z < o 



E,(Z) -E,(log(Z))- 



7/7T 



log 



1 E«(Z) 



+ log(e Ei 



TT + GiiZ) 
(log(Z))\ 



= 



?/7T 



^ + G;(Z) 



Rearranging terms, we get, 

'n + Gt(Z) 



log 



log(E,(Z)e 



-E,(log(X)) 



//7T 



7T + Gj(Z) 



= 



Denote k := log (Ei(Z)e- E '^ z ^). Then, we have, 
n + Gi(Z)\ 7/7r 



l0g V 
From ( fTSI ). we have 

GjOg) 

7T+Q(Z) 1?{Z <l) I 

which, with the previous equation, yields 



Gi(Z) E,(Z) 







log 



I 



HZ) 



+ h - f] 1 



z 







Recall that Z is actually Z (7r) . We now find that tt appears 
in the equation only as l(ir). Hence we can view this as 
an equation in the variable Z(= j). Rearranging terms, we 
get 



'E,(Z)\ E t (Z) 
- log —r- + V- 



l 



I 



-{h-n) 



Exponentiating both sides, and substituting back for bi, 
yields 

On cancelling Ej(Z), and transposing terms, we next 
obtain 

(M£l_!) = E;(log(f )} 



or. 



-^(¥)) =fi -E,Mf)) 



Taking log on both sides, we have, 



r]Ei 



Z-l 



i ;= Ei ( l0 ^(f 

In terms of G;(-), this is equivalent to 

*(¥)-«Nt 

which is the desired result. ■ 
We next address the question of a unique positive 
solution of ( fTol l. The following lemma guarantees the 
existence of a unique positive solution, when /(•), the 
density of satisfies a certain monotonicity condition. 
Lemma A.2: ( fTol l has at most one positive solution if for 

any < l\ < 1%, \ v [ 2 1 is a strictly monotone decreasing 

function of y. 

Proof: Expanding G;(-), fWi becomes, 



') I IT -1 ) 9{z)dz - / log (j) g(z)dz = 



Rewriting the equation in terms of /(■), we have, 

/Mi- 1 )-*®)?'®*- 

Using a substitution y = f in the above equation, we get, 

J Q (io g (») - V (y i))^/ (~) rfy - o (19) 

Define c(y) := (log(?/)-r?(y-l))^ and 6;(y) := / (ij. 
We are now interested in a positive Z that solves 

/ c(y)bi(y)dy = 
Jo 

Observe that lim^o c(y) = — oo and c(l) = 0. Further, 
there exists a unique y' such that c(y) < for all < y < 
y' and c(y) > for all y' < y < 1. Since i>;(y) > for 
all y and Z, we have c(y)bi(y) < for all < y < y' and 
c(y)bi(y) > for all y' < y < 1. In particular, 

c{y)bi(y)dy < 



c(y)bi(y)dy > 



y' 

Consider ii, Z2 such that < Zi < Z2. By hypothesis, 
^ 2 ^ is a strictly monotone decreasing function of y. 

Hence, ; :, 2 ) [ is also a strictly monotone decreasing 
function of y. We then have, 

So Hv)hMdy Jo Hv)\ b i^K(y)dy ^ W) 
J y ' \c(y)\b h (y)d y So' \<v)hMdy > Mt/0 ' 

And, 

S„,c(y)b l2 (y)dy Sy><y) b T$$bldy)dy ^ bj Jg) 

" My') 



< 



S v > c (y) b h (y)dy Sy> c iy) b h (y) d v 
So\<y)\biMdy Sy>c(y)t>iMdy 



>y 
Hence 



St \c(y)\b h (y)dy > Sy> c (y) b h(y)dy 

Interchanging terms, we get, 

S? \c(y)\b l2 (y)dy . So' Hy)\biAy)dy 
S y > c (y) b h (y)dy S y > c iv) b h iv)dy 

i.e., the ratio of the negative area of the integral to the 
positive area of the integral is a strictly monotonic function 
of Z. Hence, as Z increases, the integral ( fT9l ) can cross at 
most once, or, there exists at most one (non-trivial) solution 
for (29). ■ 



B. Proof of Theorem \5.1\ 

In this section, we will use the variables and equations 
from the discussion in Appendix lAl 

Lemma B.3: dxY y%f) is continuously differentiable 
with respect to d. 

Proof: Recall that tt := -p r . T [tt) and A(7r) (equa- 
tions (V2\ and ( fT3b ) are continuous function of tt, and tt 
itself is a continous function of d. Hence, from dlOl l. we 
see that dxT ("pr) i sa continously differentiable function 
ofd. ■ 

Lemma B.4: If H (or equivalently X :— -^f-) has a 

finite mean, then lim^o d x T (§r) = 0. 
Proof: Consider dTSt 



(I — z)g(z)dz = tt 



where I is in fact 1(tt). Talking I outside the integral, we 
get, 



I (l - g(z)dz = tt 



Rewriting the integral as an expectation, we have 



ZE 2 (l-f) + =7ror E,(l-f) 
tonce Convergence Theorem, we get, 

lim E z ( 1 - f 



Using Mono- 



or, 



lim — 

l — >oo / 



oo (d -> 0). 



(20) 



From (fT5l l. we see that, I — > oo as tt 
Hence, we have, 

HmM = 1 

7T — >oo TT 

Now, consider the following limit, lim^o d x r(7r(d)), or 
equivalently, lim^^oo tt~ v T(tt). We know that, 

7r~^r(7r) > 



From ( TT4b . we have, 

7r _ 'j'r(7r) = 7r"E 2 



log 



1{tt) 



TT 1 E, 10 



+ log(7r) 



Expanding the term inside the expectation, we have, 

Using the inequality log (-) < - (for z > 0) in the above 
inequality, we get, 

< tt-^E z + log ( l -^pj + log(7r; 



E z (-^) < oo (follows from the definition Z := ^ and 
the hypothesis on EX), -q > and from d20"l ). we have the 
right hand side of the above expression — > as tt — > oo, 
which implies that linijr^oo Tr^T^) = 0, or 



lim d x r(7r(d)) = 



Lemma B.5: Let w > 2, (1) be continu- 

ously differentiable and lrnx^o tt/ ( ) = 0. Then 

ro( dxr (^)) <0asd^oo.* 

Proof: From dTOb and the discussion in the proof of 
Lemma |A. 11 we have, 

d 



dd 



(d r(7r(d))) = r(7r(d)) - r77r(d)A(7r(d)) 



-j( (- (7 - 1 ) -"(f)) 

where k > 0. Using a substitution 77 = f, we get, 

i(dr(7r(d))) 



1 



dy (21) 



Define := 77 (j/ — 1) — log(y). For 77 > 1, there exists 
a y' (depending on 77) such that 6(77) > for < y < y' 
and b(y) < for y' < y < 1, also 6(1) = 0. Then, in ( Eli 
we see that, 



( v (y-l)-log(y))^f(J-)dy>0 
J (77(2/ - 1) - Iog(y))^/ (i) dy < 



Further, 



(77(2/ - 1) - log(y))dy = 1 - - 
z 



For 77 > 2, the integral J Q b{y)dy is non-positive. 

Let (7(77) := p-/ f^Y Then is continuously dif- 
ferentiable function and limj,_>o 9(77) = (by hypothesis). 
Define 770 as 

770 := sup{7/ : g(z) = 0, < z < y} 

If 770 > 0, then, we see that for I sufficiently small, 

jf" (77(77 - 1) - log(7/))-^/ (i) dy = 

This is because for sufficiently small I, ^r/ fp) — for 
< 77 < y 1 . Hence, lim^ ^(dxT (%)) < 0. 



If yo = 0, we then have g'(y) > in a small 
neighbourhood of (since g is continuously differentiable 
by hypothesis). Hence, the function g(y) is a monotonic 
increasing function in an e neighbourhood of 0, i.e., 
5(0) < g(y) < g{y') < g(e) for all < y < y' < e. 
Hence for all sufficiently small I, is a monotone 

increasing function of y in [0, 1]. Hence, in (fSTJ, we have, 



( v (y-l)-\og(y))-^f(^ )<!,, i 



y 2 l 2 yyl^ 
j\{y-l)^o g{ y))±- 2 f{^ d y 



^ 2f {^C {7]{y - i] ^ og{y))dy 



-)f(- 

y'lj J \y'l 



1 - 



V 



The final expression is non-positive for rj > 2. Thus, 
|j(dxr(|v)) <0asd^oo. ■ 

Lemma B.6: Let r\ > 2 and 4j f (j) be continuously 
differentiable. If for large x, P(X > x) = 0(4j) (or 

equivalently for H = ^-), then lim^oo d x T {^r^j = 
Proof: Let P(X > x) = O(^) for large x. i.e., 

r f(x)dx=o (\ 

J X \ x / 

Using a substitution z = -, we have, 
1 „ fl 



(22) 



Since g(z) > and continuous (by hypothesis), we have, 
<?(0) = 0. Suppose not, then, we have g(z) > e for all 
< z < 6 for some 8. Then, 



Define g{z) := 4,/ (i). Then, 

3 (z)dz = 0(z 2 ) 



g{z)dz > ez 



for all z < <5, which is a contradiction to ( 1221) . Hence 
lim 2 ^o50) = or lim 2 ^ 4^/ (A) = 0. 
We know from (TTOb that 

^(d r(vr(d))) = I>(d)) - »77r(d)A(7r(d)) 



Now from Lemma IB. 51 we see that, for rj > 2, and for 

d — > oo, 

r(?r(d)) - ry7r(d)A(7r(d)) < 
In other words, 

r(7r(d)) < r;7r(d)A(7r(d)) 
Multiplying by d on both the sides, we have, 

P' 

d T(n(d)) < W (d)X(n(d))d = r,-^\{i,{d)) (23) 

Since ^ (d T f^J")! — as d — > oo, the function 
d r(7r(d)) is monotonic decreasing for d — > oo. Also 
d r(7r(d)) > 0. Suppose that, lim^oo d T(ic(d)) ^ 0, 
it implies that lim^oo d T(w(d)) > e > 0, which, using 
( l23l . implies that > e or as d — ► cx) 

A(tt(<0) > ed^ 1 

From ( fT3l l. we have, 

ignoring the negative term, we have, 

1 f°° Pt 
f(x)dx > — — 



(24) 



d'' 



or, 



poo p * 

J f{x)dx > -±\ 



Substituting from (l24l . we have, 



P 



f(x)dx > ^-tdP- 1 = Pte\ 
d^ d 



(25) 



But we have 



f /(*)<** - P(AC > A) = O (jp) <o[^-_ 

(26) 

where the last inequality follows from (|24t . For rj > 2, 
( f25l ) and ( f26l i yields a contradiction. Hence, lim^oo d x 

r(§)=G. " 



C. Discrete Fading States 

The optimization problem for the discrete fading 
state case, simplifies to 



max ah In I 1 



hen 



P(h) \ 
dn J 



subject to a,hP(h) < P t ' 



(27) 



hen 



a hi = cti, Xi = and ^ 



For notational convenience, let us index the set of fading 
states, Ti, in descending order by the index i, 1 < i < \H\, 
i.e., hi > hi > h,3 > ■ ■ ■ . Further, denote 

-,and & = 

a* 

Also, denote 

We will later recall that, for each power constraint P t ', II 
is a function of d. Using this new notation and change of 
variables, we obtain the problem 



max ^2 a i m (1 + x i£i) 

i 

subject to 2_. a i£i — n 



(28) 



We have the maximisation of a concave mapping from 
Rl^l to R subject to a linear constraint. The KKT condi- 
tions are necessary and sufficient, and the following "water 
pouring" form of the optimal solution is well known. There 
exists A (IT) > 0, such that, for 1 < i < \H\, 

+ 



/ i 



i 



VA(n) x 

with A(II) being given by 

/ 1 1 



E 



{* : x?fr)>i} 
Defining, for 1 < k < \H 



VA(n) 



n 



a k , and a k = ^ 



and n = 0, II| W | = oo, the Lagrange multiplier, A(I1), is 
given by 

-l 



A(n) = ( _L ( ttfc + n; 

\Pk 



(29) 



for II fe _i < n < n fc when 1 < k < \H\ - 1, and for 
II|7^|_i < II < oo when k = \H\. Here the break-points 
Ilfc, 1 < k < \H\— 1, are obtained by equating the values of 



A(II) on either sides of the break-points, and are expressed 



' a k + i _ ak 

Pfc+l Pk 
1 1_ 

Pk Pfc+l 



The denominator of this expression is clearly > 0, and a 
little algebra shows that, since x k +i > Xi, 1 < i < k, the 
numerator is also > 0. 

For each II, let us denote the optimal value of the 
problem defined by d28j by r(n). We infer that 



an 



A(n) 



Now, fixing the power constraint P t , and reintroducing the 
dependence on d, we recall that 11(d) = and hence 
conclude that 

Define do = oo, di w i = 0, and, for 1 < k < \ fi\ — 1, define 



Pk Pfc+i 

"fc + i _ otk 
Pfc+l Pk 



Note that = rfi^i < di^i-i < ■ • • < g?2 < d\ < d = oo. 
Now, substituting for A(II(d)) from d29l and integrating, 
yields the following result 

Theorem C.l: For given P t ', the optimal value T(d) of 
the problem defined by (l27b has the following characteri- 
sation. 

1) The derivative of T(d) w.r.t. d is given by 



dr 

dd 



VPkPt 



1 



d \ a k di + P t 



(30) 



for d k < d < dk-i when 1 < k < \H\ — 1, and for 
< d < d\n\-i when k = \H\. 

2) ^ is a negative, continuous and increasing function 
of d. In particular Y(d) is a decreasing, and convex 
function of d. 

3) The function T(d) is given by 



T(d) = p k In ( a k + j j k 



(31) 



for d k < d < e?fc_i when 1 < k < \Tl\ — 1, and for 
< d < (ii^i-i when k = \H\, with the constants 
of integration j k being given as follows. 



7i 



1 

Oil 



Xx 

ffll 




Fig. 3. The stationary points of V(d)d lie among the intersections of the 
curve ye~ nv and lines e bk ~ 1, 1 < k < \H\, in the interval < y < 1. 
Here the plot is drawn for r\ = 3. 



and, for 2 < k < H, j k is obtained recursively as 



Ik 



Oik-l 



Pt' 



7fc-i 



' Pk-i 



Proof: OTb is obtained by integrating the derivative 
in d30l > over each segment of its definition. The integration 
constants j k are obtained by equating T(d) on either sides 
of the break-points of the argument d. ■ 
1 ) Optimisation over d: Using Theorem IC.ll we con- 
clude that we need to look at the stationary points of T(d)d. 
To this end, consider the solutions of 

T(d)+dT'(d) =0 

Reintroducing the variable II = -p r , and canceling p k , we 
need the solutions of 



, ,, n\ rjU 
In H a k jk — pf 







for n k -i < n < n fc when 1 < k < \H\ - 1, and for 
_ 1 < II < 00 when k = \fi\, with the break-points 
Ilfc, 1 < k < \7i\, as given earlier. Let us write a = 
1— — ^ , define b k = lnafc7fe (observe that b\ = 0), and, 

for given k, use the new variable 

1 

V 



1 



Note that, for < II < 00, 1 > y > 0. Define S k = — \r-. 
Then we seek the solutions of 

In - + b k - n (1 - y) = 

y 

for S k < y < 5 k -i, for each k,l < k < \H\; note that 
Sq = 1, and = 0. The equations can be written more 
simply as 

e bk-v =ye -vv } 

and are depicted in Figure [3] At this point we can conclude 
the following 



Theorem C.2: There are at most 2\H\ — 1 stationary 
points of T(d)d in < d < oo. 

Proof: The result follows from the arguments just 
before the theorem statement, since each line e^ bk ~' n ' , for 
2 < k < \fi\, has at most two intersections with ye^ r,v , 
in < y < 1, and e~ r ' has only one such intersection. ■ 

D. Fixed Transmission Time vs Fixed Packet Size 

In this section, we will formally establish that fixed 
transmission time schemes are more throughput efficient 
compared to fixed packet size schemes, for a given average 
power constraint. We will prove this result in a general 
framework, without explicitly modelling the underlying 
MAC, the power control schemes used or the channel 
fading distribution. 

Data Transmission Model: In a fixed transmission time 
scheme, all data transmissions (with positive rate) are of 
a fixed amount of time T, independent of the channel 
state h and the power used. Earlier, in our work (see 
Section IIII-Bb . we assumed that, when the channel fade 
is poor (and hence P(h) = 0), the channel is left idle 
for the next T seconds. Further, the optimal power control 
policy for such a system was found to be a non-randomized 
policy, where a node transmits with constant power P(h) 
every time the channel is in state h (see Section IV- Al l. 
Here, we will allow the possibility of the channel being 
relinquished when bad with a fixed time overhead < T. We 
consider a spatio-temporal fading process with successive 
transmitter-receiver pairs being selected by a distributed 
multiaccess contention mechanism. Hence, relinquishing 
the channel might improve throughput, as successive fade 
levels might have little correlation. The optimal policy for 
such a MAC could be a randomized policy. Hence, we 
will allow a randomized power control, i.e., for a channel 
state h, the transmitter chooses a power Ph according to 
some distribution. In a fixed packet size scheme, all data 
transmissions (with positive rate) carry a fixed amount of 
data L independent of the channel state h and the power 
control used. Here as well, we will allow the possibility 
of a randomized power control and the possibility of 
relinquishing the channel with a fixed time overhead (when 
the channel fade is poor). 

Optimality Criterion: The throughput optimality of a 
data transmission scheme is established either by compar- 
ing the energy required to send a certain amount of bits in 
a given time or by comparing the amount of bits sent with 
a given amount of energy in a given time. (We will discuss 
more about this optimality criterion in Remark TD-U . We 
study a data transmission scheme by considering two data 
transmissions of positive rates, in some arbitrary channel 
states with gains hi and h 2 and with applied powers 



Phi and Ph 2 - We do not make any assumption on the 
probabilities of hi and 1%2, and about the power control 
policy which yields the powers Ph ± and Pf l2 . 

For a given power control scheme (h, P^), we will 
then assume that the transmission rate given by Shannon's 
formula is achieved over the transmission burst; i.e., the 
transmission rate is given by 

C h = Wlog(l + hP h ) 



We have absorbed the factor 



in to the term h (since 



d is fixed in this discussion). Hence, the time durations 
taken to transmit the L bits during the channel states hi 
and h 2 (with the powers P^ and P} 12 ) are given by Ti ll := 



andTL := 



Wlog(l+hiP hl ) ' Wlog(l+h 2 P h2 ) 

time occupied by these two transmissions is 



Then, the total 



L 



W\og(l + hiP hl ) W\og(l + h 2 P h2 ) 
spending an amount of energy equal to 

LP hl LP h2 



(32) 



(33) 



Wlog(l + hiP hl ) Wlog(l + h 2 P h2 ) 

Define Lp := 2 x L as the amount of bits sent in time Tp 
using an energy Ep in channel states hi and h 2 . 

Lemma D.7: Let hi > h 2 . For a fixed packet size 
scheme, if P^ and Ph 2 are applied powers during channel 
states hi and h 2 , then having hiP} l:L > h 2 Pf l2 is through- 
put optimal. 

Proof: Suppose that hiPf ll < h 2 Ph 2 . Then, 

\og(l + hiP hl ) <\og{\ + h 2 P h2 ) 

Find power controls Ph x and Ph 2 such that 

\og{l + hiP hl ) = log(l + h 2 P h2 ) (34) 
log(l + h 2 P h2 ) = \og(l + hiP hl ) (35) 

or, equivalently, 

h x P hl = h 2 P h2 (36) 
h 2 P h2 = hiP hl (37) 

With the power control scheme (hi, P^), (h 2 , Ph 2 ), the 
total time occupied in the transmissions of 2 x L bits during 
the channel states hi and h 2 is, 

L L 



Tp = 



W\og(l + hxP hl ) W\og{\ + h 2 P h2 ) 



(from OH and d35ll). Now, consider the energy spent to 
transmit these 2 x L bits, i.e., 



E p = 



LP h 



LP h 



Substituting for Pf ll and i\ 2 from d36T ) and d37l ), we have, 
1 Lh 2 P h2 1 LhiP hl 

p hiW\og(l + h 2 P h2 ) h 2 Wlog{l + hiP hl ) 
Rearranging the terms, we have, 

1 LhiP hl 1 Lh 2 P h2 



E s 



h 2 W log(l + hiP hl ) hi W log(l + /i 2 Ph 2 
1 LhiP hl 1 L/i 2 P^ 2 



< 



< 



Wlog(l + h 1 P hl ) Wlog{l + h 2 P h2 ) 



hi W l0g(l + /liP /ll ) /l 2 l0g(l + h 2 P h2 

LPh! LP h2 

W\og{\ + hiP hl ) W\og{l + h 2 P h2 ) 
= E P 

where the inequality follows from the fact that 

LhiP hl (1_ 1_ 
W\og(l + hiP hl )\h 2 hi 

Lh 2 P h2 (\ 1_ 

W\og{l + h 2 P h2 ) \h 2 hi 

since hi > h 2 and hiP^ < h 2 Ph 2 (by assumption) and 
the fact that log (^ +x ^ is strictly monotone increasing. 

It follows that an optimal power control must have 
hiP hl >h 2 P h2 . m 
Remark: From Lemma lb. 71 we see that, when hi > h 2 , 
C hl := W\og{l + hiP hl ) > H^log(l + h 2 P h2 ) C h . 2 , 
or equivalently, Th x < 7) l2 . ■ 
We will now provide a comparison of the fixed packet 
scheme with a fixed transmission time scheme and show 
the optimality of the fixed transmission time schemes. The 
comparison is done under the following assumption. 
• The channel has the same marginal fading distribu- 
tion, whenever sampled by a transmitter, for either 
schemes. This is a reasonable assumption as we con- 
sider spatio-temporal fading, with successive trans- 
missions from possibly different source-destination 
pairs chosen by the distributed multiaccess contention 
scheme. 

For the fixed packet size scheme, hp := 2 x L bits were 
transmitted in Tp(= Th t + TX 2 ) time (see ( [32l ) with an 
amount of energy equal to Ep (see (|33l), in two channel 
samples hi and h 2 . A reasonable comparison would be to 
find the throughput of a fixed transmission time scheme 
for a total duration of Tp seconds involving two data 
transmissions with channel samples hi and h 2 of equal 
duration T = ^ and a total energy of Ep. We will assume 
that P/ ll and P} 12 , the power used for the fixed packet size 
scheme are such that Th 1 < Th 2 (see Lemma II). 7I >. Hence, 
we have T\ 1 < T < T/ l2 , or, the fixed transmission time 
scheme spends relatively more time on a better channel. 
Clearly, its throughput is better than the fixed packet size 
scheme for the same energy constraint, as seen below. 



Let Pt h and Pt h be the optimal power control for the 
fixed transmission time strategy such that 



E T := Pt h T + P th T = P hl T hl + P h2 T h2 = E P 
We have, 

L P = 2L = T hl W \og{\ + hM + T h2 W \og(l + h 2 P h2 ) 
Expanding the left hand side, we have, 

2L = T hl W\og(l + h 1 P hl ) + (T h2 ~T)W\og(l + h 2 P h2 
+ TW\og(l + h 2 P h2 ) 

Using hi > h 2 , we get, 

2L < T hl log(l + hiP hl ) + (T h2 - T) log(l + hxPhJ 
+ Tlog(l + h 2 Ph 2 ) 
< TlogCl + ZuP^J+TlogCl + ^J 
=: Lt 

where the last inequality follows from the fact that 
(hi,P thi ) and (h2,Pt h2 ) is the optimal power control 
scheme for the fixed transmission time scheme with time 
T P (= 2 x T) and energy E T {= E P ). 

Remarks D.l: For L(t) defined as the amount of bits 
sent up to time t, and E(t) defined as the total energy 
spent up to time t, the average throughput (6) and the 
average power (P) of the system are, in general, defined 
as 

6:=liminf^ 

t— too t 

E(t) 

P := lim sup 

t — >oo i 

Under additional assumptions on the fading process and 
the power control scheme used, the expressions are sim- 
plified as an ensemble average (for example, see ([U and 
d2J for a fixed transmission time scheme). In this section, 
the optimality of the schemes have been shown directly, by 
comparing the amount of bits transmitted for a particular 
sample of channel for a given amount of time and energy, 
or by comparing the amount of energy used to transmit a 
given amount of bits for a particular sample of channel in a 
given amount of time. For example, the argument provided 
here directly translates to an argument with the ensemble 
average for the discrete fading case. This approach is not 
only straightforward, but also is very general. ■ 



